Analysis of D^^ipn^ beyond naive factorization 

Haijun Gong ^ Junfeng Sun ^ Dongsheng Du "'^ 

a. CCAST (World Laboratory), P.O.Box 8130, Beijing 100080, China 
h. Institute of High Energy Physics, Chinese Academy of Sciences, 
P.O.Box 918(4), Beijing 100039, China ^ 

February 1, 2008 



Abstract 

We analyze the decay Ds^(f)iT with QCD factorization in the heavy quark hmit. 
The nonfactorizable contributions, including hard spectator contribution are dis- 
cussed and numerical results are presented. Our predictions on the branching ra- 
tio of the decay are in agreement with the experiment. We also use a pure phe- 
nomenological method to estimate the branching ratio for Ds^(j)TT with the existed 
D^^K*Tr data. 
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1 Introduction 



Both CLEO [|T| and BES [§] have reported their direct model-independent measurements 
for the Dg— i>07r branching fraction: 



Br(D ^a,^.ii^-59±0.77±0A8)xl0-' CLEO 
br{V, ^ 07rj - I (3 9+5.1+1.8) ^ g^g^ 



The average branching ratio of D^— >07r is (3.6 ± O.9)xlO~2[0]. 

The precise estimation of the branching ratio for the decay Ds^4'Tt is very important. 
First, it is difficuh to measure the absolute branching ratio of Ds^(j)TT because we do 
not know the fraction of D^D~ pairs production in e^e~ annihilation in comparison 
with DD pairs ( BES used e^e~ — > D^D~ to obtain the first direct model-independent 
measurement of the Ds—>-(j)7r branching fraction, however, with only two "double-tagged" 
events ). But we need to know the branching ratio for the study of B decays such as 
B^DgX etc. Moreover, most of the measurements of the meson branching fractions 
are normalized to the clean Ds^4>'n' channel. Second, theoretically, the decay of Ds^4>7r 
is dominated by spectator diagram with external emission of pion. This is easier to handle 
compared with other exclusive non-leptonic decay channels. 

Previous calculations for the branching ratio Br{Ds^(f)TT) are based on the naive fac- 
torization approach which is proposed by Bauer, Stech and Wirbel(BSW) |Q. But in 
BSW approach, non-factorizable effects can not be calculated, they have to be parame- 
terized by an effective color number A*"^-^-^ which is treated as a free parameter. Moreover, 
results obtained with BSW approach still depend on renormalization scale and scheme. 
The authors in [^] examine the Ds—>-(f)'K amplitude through a constituent quark-meson 
model. With this model, the calculated decay width r{Ds—>-(f)TT) is larger than the exper- 
imental data. Paver and RiazuddinQ studied Ds^cpn in a valence quark triangle model, 
incorporating chiral symmetries, the result is compatible with the experimental data. In 
1]^, ^, the authors considered the contribution from the color octet: {(j)7i~^\H^\D'^), where, 
H^^ = lEa {uX^c){s\''d). But they all introduced some new parameters, so they brought 
new theoretical uncertainties. 

In the past years, Beneke, Buchalla, Neubert and Sachrajda developed QCD factoriza- 
tion (QCDF) approach |^ to calculate the hadronic matrix elements of B decays in the 



heavy quark limit. It has been used for many B decays modes [^, ^ with interesting 
results. In the present paper, we will follow this method to calculate the branching ratio 
for Ds—>-(f)TT. In the heavy quark limit itLc'^^Aqcd, non-factorizable contributions are con- 
sidered from the first principle. In Ds—>-(j)'K decay, the hadronic matrix elements can be 
represented as: 

(7r0|g,(/^)|D,) = (7r|Ji|O)(0|J2|D.)-[l + E^„< + O(AQC7D/me)]. (1) 

The naive factorization corresponds to neglecting the 0{as) corrections and the power 
corrections in Aqcd/^c- Although rric is not as large as m^, we still hope that the QCD 
factorization approach in the heavy quark limit can also give a reasonable description of 
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Figure 1: Dependence of the Wilson Coefficients Ci, C2 on the renormahzation scale /i at leading 
order (a) and next-to- leading order (b). 



Ds meson hadronic decays. With this method, we analyze Ds^4>t^ decay and compare 
the results with those obtained with naive factorization. Finally we use existed data of 
D^^K*~TT^ to estimate the branching ratio Br{Ds^4'7!') in a model independent way. 



2 Ds^cjyn in QCD Factorization 

The low energy effective Hamiltonian for Ds—>-(j)7T can be expressed as follows : 

^e// = ^v:^vuCiif^)Qiif^)+C2if^)Q2if^)]. (2) 

The four-quark local operators Qi^2 are 

Ql = {SaCi3)v-A{Ul3da)v-A ^g-j 

where a, j3 are the color indices of SU (3)c- Wilson coefficients Ci(/i) are universal, process- 
independent and calculable with the renormalizaion group improved perturbative theory, 
their /x-dependence are expected to be cancelled by the hadronic matrix elements. The 
leading order (LO) and next-to- leading order (NLO) corrections to Cj(/i) have been pre- 
sented in [0. In the naive dimensional regularization (NDR) scheme, we give the numer- 
ical values for Ci{fi) at three renormahzation scales in Tab.l. In Fig.l, we also display 
the dependence of Ci{i = 1, 2) on /i in the LO and NLO approximation. We will take the 
values of Ci{i = 1, 2) at NLO for our forthcoming calculations. 
The decay constant and form factors are defined by : 

(7r+|(nrf)f^_^)|0) = -z/.p^, (4) 
msc),iv-A)\Dt) = -^{{mn.+m^)e;A?^%') '-^{pd. + P^),A^'%') 



Table 1: Wilson coefficients in NDR scheme. The input parameters in numerical calcu- 
lations are fixed: as{mz) = 0.1185, aemijnw) = 1/128, my/ = 80.42 GeV, niz = 91.188 
GeV, mt = 168.2 GeV, rUc = 1.45 GeV . 



/i 




C^2(/i) 


LO 


NLO 


LO 


NLO 


= 


IGeV 


-0.650 


-0.500 


1.356 


1.272 


= 


rric 


-0.520 


-0.390 


1.268 


1.200 


= 


2GeV 


-0.435 


-0.314 


1.214 


1.153 



+ ^- (5) 

where = {pos -P4>)t^ and 

<^^(0) = A?=^(0), (6) 
2mX^'^(0) = (m,5.+m^)Af»^(0)-(mB, -m^)Af='^(0). (7) 

The relations (6) - ensure that there is no kinematical singularity in the matrix element 
at g2 = 0. 

Under naive factorization, using Eq. (4) - (|^, the decay amplitude of Dg-^tpTT reads 



A{Ds^<Pti) = V2GFV:,Vudhm^A^^'^{ml){e* ■ poj ■ a^, (8) 

where 02 = C2 H ^t/Ci; ^^c^"^ number of colors. The form factor Aq'''^ is defined 

by d^). From Eq. (|^) we can see that the amplitude depends on the renormalization scale 
/i, because the Wilson coefficients Ci{n), C2(/i), and hence oi, a2 depend on /x, whereas 
the decay constant and form factor are independent of fi. So the amplitude A{Ds^(f)n) 
is /i-dependent. On the other hand, it does not consider the nonfactorizable effects. If 
we calculate it with QCD factorization, take all the high order corrections into account, 
Qi and the amplitude A{Ds^(l)n) will be fi independent. In our paper, we calculate it 
only to the order of as, so and the amplitude ^(Ds— >-07r) still depend on /i, but the 
dependence is less sensitive to fi. With these preliminaries, we now analyze Ds^4>7i' with 
QCD factorization. 

In Ds-^(f)7T decay, the emitted meson vr is light, the hadronic matrix elements can be 
written as: 

{7r<P\Qi{fi)\Ds) = A^^^{{))CdxTl{x)^^{x) 

J 

+ fdi dx dy Tl\i,x,y)^DXO^A^)^M- (9) 
J 
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^^"''^(0) denotes the Dg^cj) transition form factor, $d(0) ^-wi^) and label light- 

cone distribution amplitudes (LCDAs) of Dg, n and (p meson respectively. Tj ' de- 
note hard-scattering kernels which are calculable in perturbative theory. Neglecting 
the 0{Aqcd/^c) corrections, T/'^^ are hard gluon exchange dominant. Other non- 
perturbative contributions are contained in the LCDAs of mesons or the form factor. 
The second term in Eq. represents the hard spectator contribution. 

We next proceed to calculate the nonfactorizable effects in the D'^^(j)7c~^ with QCDF 
approach. Then in heavy quark limit, for simplicity, we will neglect the masses of light 
quarks and vr. We consider the vertex corrections and hard spectator interactions depicted 
in Fig. 2. The technique is similiar to that of the B^nir/K mode, readers can be referred 
to for details. As in 0, we obtain the QCD coefficients aj(i= 1, 2 ) at NLO in NDR 
scheme. Then the coefficients are given as 

Here = 3 (/ = 4) is the number of colors (flavors), and Cp = ^2n^ is the factor of 
color. We define the symbols in the above expressions as the same as Beneke's, which are 

F = -18 -121n-^ + + (11) 

rric 

fi = [dxg{x)^^{x), (12) 
Jo 

with the hard-scattering function 

1 - 2x , 

g{x) = 3 ma; — Sin. 

1 — X 

The hard spectator scattering contribution is given by 



- wJirr- /^«%^ /^-^ f'y'-^- (13) 

N ^(0)m|)^ ^0 ^ Jo X Jo y 



where (/dJ is the (Dg) meson decay constant, m^^ the mass of Dg meson, Aq'''^{0) 
the Dg ^ (j) transition form factor at zero momentum transfer, and ^ the light-cone 
momentum fraction of the spectator quark in the Dg meson, /// depends on the wave 
function through the integral 

f di^oAO/i = ^Dj\D. (14) 

Jo 

This introduces a new hadronic parameter Ad, Ad is of order Aqcd, we take A^ = 335 
MeV here. 



5 



XX XX 



(a) (b) (c) 




(d) (e) (f) 

Figure 2: Order of as corrections to the hard scattering kernels T/ and T^^ . The two hnes 
directed upwards represent the two quarks that make up vr. These diagrams are cahed vertex 
corrections for Fig.(a)-(d) and hard spectator diagrams for Fig.(e)-(f) respectively. 



From the expression (10) of the QCD coefficients ai{i = 1, 2), with the renormahzation 
group equation for Wilson coefficients Ci{fi) at leading order logarithm approximation 

dCiifi) _ as ^Tr~< ^.A (^15) 



din fi 477 

where 7 is the anomalous dimension matrix, we find = (i = 1,2) at the order of 
as, this makes the /i-dependence of the decay amplitude calculated with QCDF approach 
less sensitive than that calculated with naive factorization. This point can also be seen 
roughly from the data in Tab. 2 and Fig.3 - 4. But there are still uncertainties in the 
calculation, such as the form of wave functions and unknown form factor Aq"'^. 

Notice that in the decay Ds —>■ (pTc, using the isospin analyses |1^, |T^, we find that the 
final state involves only a single isospin, so there is no interference effects from the final 
state interactions (FSI) when we calculate the branching ratio of Ds — > 0vr. 

In the Ds rest frame, the two body decay width is 



1 



T{Ds^<Pn) = —\A{Ds^<PTr)\ 



\P\ 



m 



2 ' 



(16) 



where 



\P\ 



m 



Ds 



2m 



Ds 



is the magnitude of the momentum of (p nieson. With the approximation m^/m|,^ ^ 0, 
the decay width is given by 



F {Dt ^ 07r^ 



327r 



iK.riK.ri«2p ( 

X^Ds 



.rriDs 



(<=^(0))^ (17) 



The corresponding branching ratio is given by 

T{Ds ^ 07r) 



Br{Ds ^ 07r) 



F 



total 



total 



1 

'TDs 
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In our numerical calculations, we will take the following values for the relevant input 
parameters [|: = = 0.975, U = 131 MeV, = 233 MeV, m, = 1.45 GeV, fD,= 
280 ± 19 ± 28 ± 34 MeV. As for the form factor A^''^(0), for lack of experimental data, 
we use the value taken from the Ref. 0] Aq"'^{0) = 0.70. For mass of the mesons, we use 
mD = 1968.6 ± 0.6 MeV, m^= 1019.417 ± 0.014 MeV. If not stated otherwise, we shall 
use the central values as the default values in our later calculations. 

For distribution amplitude of vr, two kinds of the wave functions are used, one is the 
asymptonic form |^ ^n{x) = 6x(l — x), the other is delta-function ^t^^x) = — |). In 
Tab. 2 we list the values of ai,a2 and branching ratio(i3r) at fi= 1 GeV, rric, and 2 GeV 
with different wave functions of vr. The numerical results which are calculated with BSW 
approach ( Where we take N^-^-l" = oo because the experimental data of MARK III for 
charm decays do not show color suppression ) are also listed for comparison. 



Table 2: The values of and Br at jj,= 1 GeV, rric and 2 GeV calculated with QCDF and 
BSW approach {N'^f^ = oo ). For QCDF, we calculate the spectator contribution with 
three forms: /// = and two different wave functions of vr. In the QCDF columms, the 
values in the parentheses are those with ^Tr{x) = 6x{l — x), the values in the brackets are 
those with $7r(a;) = S{x — |). 



/i 


ai 


a2 


Br % 


BSW 


QCDF 


BSW 


QCDF 


BSW 


QCDF 


/i= 1 GeV 


-0.500 


-0.396 - 0.215 i 
( -0.071 - 0.215 i ) 
[-0.168- 0.215 i ] 


1.272 


1.231 + 0.084 i 
( 1.103 + 0.084 i ) 
[ 1.141 + 0.084 i ] 


3.75 


3.53 
( 2.84 ) 
[ 3.03 ] 




-0.390 


-0.284 - 0.150 i 
( -0.058 - 0.150 i ) 
[-0.125 - 0.150 i ] 


1.200 


1.165 + 0.049 i 
( 1.092 + 0.049 i ) 
[ 1.112 + 0.049 i ] 


3.33 


3.15 
( 2.77 ) 
[ 2.88 ] 


IJ= 2 GeV 


-0.314 


-0.213 - 0.119 i 
( -0.032- 0.119 i ) 
[-0.086 - 0.119 i] 


1.153 


1.126 + 0.033 i 
( 1.077 + 0.033 i ) 
[ 1.091 + 0.033 i ] 


3.08 


2.94 
( 2.69 ) 
[ 2.76 ] 



It is necessary to note that the QCDF approach gives ai{i = 1,2) an imaginary part, 
which comes from the gluon exchange of the quarks u and d in n with the s quark in 
(see Fig. 2 (c)-(d)). Moreover, the imaginary parts of ai{i = 1,2) have no relations with 
From the numerical values summarized in Tab. 2, we find that the vertex correction 
in Fig.2 (a)-(d) is about 5~ 7%, the hard-spectator diagrams can reduce over 10% of the 
values obtained with BSW approach. And the coefficients ai{i = 1,2) are less sensitive 
to the choice of the wave functions. In Fig.3, we depict the dependence of ai,a2 and 
Br on scale /i when considering the vertex corrections (but neglecting the hard spectator 
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Figure 3: Dependence of ai,a2 and Br on the renormalization scale /i in BSW and QCDF(/;; 
= 0). The dotted and dashed hnes correspond to the values obtained with BSW and QCDF 
approach. 

contribution), we also show the results calculated by BSW approach for comparison. The 
horizontal solid lines in Fig.3(b) show the experimental branching ratio at la level. It is 
clear that the scale dependence of the values calculated with QCDF approach are milder 
than that with BSW approach. But the fi dependence still exists, the reason is that we 
calculate only at one-loop level, the source of fi dependence is from the high order 
effects. When considering the contributions from the high order corrections in ag or 
^qcd/ttIc, the fj, dependence of our predictions will be further reduced. 

In Fig. 4, we compare the results which are calculated with different wave functions 
of TT when considering the hard-spectator contribution in QCDF. It shows again that 
ai,a2 and Br are less sensitive to the selection of the wave function of tt, moreover, 
their yu-dependence is furthur reduced. From Fig. 3 and Fig.4, we find that the results 
obtained with QCD factorization approach fall in the la allowed region from the central 
experimental value 3.6x10"^, regardless of the seletion of the function of tt. Though 
the branching ratios with BSW approach are also within the la region, this approach 
takes N^-l'-l' = oo in order to fit the experimental data, so it is more phenomenological 
in comparison with QCDF approach. From Fig. 3 and Fig.4, we can see apparently that 
our predictions with QCDF approach are small compared with the values obtained with 
BSW approach. 

3 Direct estimation of Br^Dg^cpTv) 

Now we estimate the Br{Ds^(l)'n:) directly from the existed data of D^^K*^'k^ . As- 
suming spectator diagram dominance, Df^(j)7r~^ can go through quark decay diagram 
depicted in Fig. 5(a), the decay width is (16). Using the experimental data listed in Sec. 2, 
we get 0.720. 

Consider the decay D^^K*~7i which proceeds dominantly through diagram Fig.5(b), 
obviously, in Fig. 5, diagram (a) and (b) are very similiar, if s in (a) is replaced by u, we 
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(b) 

Figure 4: Dependence of ai,a2 and i3r on the renormalization scale ^ in QCDF with different 
function of vr. The dotted and dashed hues correspond to the values obtained with ^ni^) = 
6x(l — x) and ^Tr{x) = 6{x — ^). 

will get (b). In addition, the particle decay width of D^^K*7c is 

r(D^ K*n) = —\A(D'' ^ K*n)\^^, (19) 

where \p/\ = 0.719. The momentum of K* in the rest frame is almost the same as that 
of in D'^^(f)TT'^. So the Lorentz contraction effects of the wave function of and K* are 
nearly the same. We know that the decay amplitudes A{Ds — > (pir) and A{D^ K*t^) are 
proportional to the wave function overlap integrates of — and — K* respectively. 
Moreover, |p|= 0.720 and \pi\ = 0.719 mean that these overlap integrates are almost the 
same under the condition of SU(3) symmetry. The SU(3) symmetry breaking effects in 
the case of D^^(j)7r^ and — > K*^^^ should be fairly small. So in approximation, we 
can have 

\A{Ds 07r)|^|A(D° ^ K*-n+)\. (20) 

Using the experimental data 0: 

r(D°) = (0.4126 ± 0.0028) xlO^^^s, 
i3r(Z}° -> K*-7f+) = (5.0 ± 0.4)%, 

r(D.) = (o.mt'oZ) X 10-^^^, 

with Eq.(18) - (20), we obtain Br{Ds (pir) ~ (5.40 ± 0.45)%, where the error comes 
from that of the data of t{D^), Br{D^ K*~-k^) and t{Ds). It is a little outside the 
one a allowed region from the central experimental value 3.6x10^^. 
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Figure 5: Diagrams for the decay Ds^(/)7r and Z)° K*7r. 



4 Conclusions 

We have analyzed the decay Dg^ip'n: with QCD factorization in the heavy quark hmit. We 
calculate the nonfactorizable contributions, including vertex correction, hard-spectator 
contribution. These nonfactorizable contributions can give over 10% corrections to naive 
factorization. Moreover, according to our calculations, the branching ratios with QCDF 
approach is not sensitive to the choice of the wave function of pion. Our predictios are 
in agreement with the present experimental data. The direct estimation of i3r(Ds— >-07r) 
from D^^K*TT data gives a bit larger result comparing with the present data. But the 
measured data on Br{Ds-^(f)TT) are still rough, we need more data for drawing our final 
conclusion. 
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